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Abstract. Despite the absence of viscous drag, the neutron superfiuid permeating 
the inner crust of a neutron star can still be strongly coupled to nuclei due to non- 
dissipative entrainment effects. Neutron superfluidity and entrainment have been sys- 
tematically studied in all regions of the inner crust of a cold non-accreting neutron star 
in the framework of the band theory of solids. It is shown that in the intermediate 
■ layers of the inner crust a large fraction of "free" neutrons are actually entrained by 

Qh the crust. The results suggest that a revision of the interpretation of many observable 

astrophysical phenomena might be necessary. 



1. Introduction 

No sooner did Bardeen, Cooper and Schrieffer (BCS) publish their theory of supercon- 
ductivity than Migdal speculated about nuclear superfluidity in the interior of neutron 
stars. This possibility was further investigated by Ginzburg and Kirzhnits in 1964. The 
discovery of the first pulsars and the observation of frequency glitches followed by 
very long relaxation times of the order of months provided very strong evidence for 
the presence of superfluids in neutron stars. In particular, pulsar glitches are generally 
believed to be related to the neutron s uperfiuid that permeates th e inner crust of neutron 
stars at low enough temperatures (see lChamel & Haenselll2008l for a review). The core 
of a neutron star may also contain different kinds of superfluids and superconductors 
but their nature and their properties still remain very poorly known. The existence of 
a neutron superfiuid in neutron-star crusts has recently found support from the obser 



vation of the initial coolin g in persistent Soft X-ray Transients (Shternin et all 12007 



iBrown & Cummirig1l2009l) . Superfluidity is also expected to play a role in other astro- 
physical phenomena lik e the quasiperiodic osc illations detected in the giant flares from 
Soft Gamma Repeaters (lAndersson et al.ll2009h . 



2. Nuclear superfluidity 

Most studies of nuclear superfluidity have focused on homogeneous neutron matter. 
As the temperature T falls below some critical temperature T c , neutrons with anti- 
aligned spins and zero orbital angular momentum form Cooper pairs like electrons in 
conventional superconductors. Because these pairs are bosons, they behave coherently 
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on a very large scale. As a consequence, the neutron condensate can flow without any 
viscosity, analogous to superfluid helium-3. 

Recent microscopic calculations using realistic nucleon-nucleon interactions and 

based on different many-body methods tend to predict t he same density dependence of 

the n eutron pairing gap, at least at not too high densities (Chan g~et al.l2004l : lGezerhs & Carlso"n 
I2010I) . In fact, the analytical expression of th e gap in the limit of very dilute neu - 
tron matter has been known for a very long time (Gor kov & Mel ik-Barkhud arovlll961 ). 
These calculations predict a maximum pairing gap A ~ 2 MeV corresponding to a 
critical temperature T c = exp(y)A/(/:B^) ~ 10 10 K where y is the Euler-Mascheroni 
constant and Icq the Bolzmann constant. 

However in spite of their names, neutron stars are not only made of neutrons. In 
particular, the neutron superfluid present in the inner crust coexists with neutron-rich 
clusters arranged in a regular Coulomb lattice and with a highly degenerate relativistic 
electron gas. Because the superfluid coherence length is much larger than the cluster 
size, the pairing phenomenon is highly non-local. This means that even if the center 
of mass of a neutron pair is located outside clusters, one of the partners may actually 
lie inside so that both bound and unbound neutrons are actually involved in the pairing 
phenomenon. As a result of these proximity effects, the properties of the neutron su- 
perfluid in the crust of a neutron star may differ substantially from those calculated in 
pure neutron matter. 

We have recently studied the inhomogeneous matter of the crust using the nuclear 
energy density functional (EDF) theory dChamel et al.ll201Qf) . This theory, which has 
been already very successfully applied to the description of medium-mass and heavy 
nuclei, provid es a consistent des cription of both the neutron superfluid and the nuclear 
crystal lattice dChamel et all 120121) . For this purpose, we have made use of accurately 
calibrated nuclear EDFs based on Skyrme effe ctive nucleon-n ucleon interactions com- 
pleme nted with a a micros copic pairing EDF (Chamel bolOh . The Brussels-Montreal 
EDFs (iGoriely et al.ll2010h yield an excellent fit to essentially all experimental atomic 
mass data (2149 nuclei) with rms deviations lying below 0.6 MeV, while reproducing 
at the same time i) a realistic neutron-matter equation of state and ii) the 1 Sq pairing 
gaps both in symmetric nuclear matter and in pure neutron matter, as obtained from 
many-body calculations including medium polarization effects. For these reasons, the 
Brussels-Montreal EDFs are well-suited for studying the pairing phenomenon in the 
inner crust of neutron stars. The neutron-star crust calculations have been carried out in 
the framework of the band theory of solids. In this way, the effects of the periodic crys- 
tal lattice on the neutron superfluid are properly taken into account. Because the pairing 
mechanism directly arises from the strong attractive nuclear force, many neutrons are 
involved and not only those lying on the Fermi surface as in conventional supercon- 
ductors. Neutrons belonging to different bands feel different pairing interactions thus 
leading to a dispersion of the pairing gap of a few hundred keV around the Fermi level. 
Moreover, pairing correlations are found to be enhanced in clusters while they are re- 
duced in the intersticial region. These proximity effects become very important as the 
temperature approaches T c , as illustrated in Fig.Q] All in all, the critical temperature is 



found to be reduced by ~ 10 - 20% as compared to pure neutron matter dChamel et al. 



2010). The effects of the crystal lattice on the superfluid dynamics are much more 



dramatic. 
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Figure 1 . Neutron pairing field in the Wigner-Seitz cell of the inner crust of a 
neutron star at the average baryon density h = 0.06 fin -3 and for T = 0.87V Solid 
line: band structure calculation. Dashed line: at each point in the cell, the pairing 
field is approximate d by its value in pure neutron matter for the corresponding local 
neutron density. See Chame l et al.1 (1201 2l) for further details. 



3. Entrainment 



Even though t he neutron superfluid can flow with out friction, it can still be entrained 
by the crust dCarter et al.ll2006l : iPethick et al.ll2010h . Unlike viscous drag, this entrain- 
ment effect is non-dissipative. It arises because neut rons can be elastical l y sca ttered 
by the nuclear latti c e. Sin ce the first experiments by Hal ban & Preiswerkl (1 19361) and 
Mitchell & Powers! (119361 ) a few years after the discovery of the neutron, the neutron 
diffraction has been routinely applied to study crystal structure. This is because neu- 
trons reflected by parallel lattice planes spaced d appart can interfere constructively 
but only for specific incident angles 9 (measured from the plane) given by Bragg's law 
2d sin 9 = NA, where N is any integer and A < 2d is the de Broglie wavelength of 
neutrons. Similarly, neutron diffraction can occur in the inner crust of a neutron star. 
Note that Bragg diffraction is a highly non-local phenomenon. Whereas neutrons used 
for diffraction experiments are produced in nuclear reactors and accelerators, the neu- 
trons in the crust of a neutron star originate from electron captures on nuclei during 
the gravitational core-collapse of the supernova. With increasing matter density, nuclei 
thus become more and more neutron rich. As the matter density reaches about 4 x 10 11 
g cm -3 , some neutrons drip out of nuclei and form a highly degenerate quantum liquid. 
A "free" neutron that is Bragg reflected cannot propagate and is therefore trapped in the 
crust: its velocity in the crust frame is zero. In other words, this neutron is entrained 
by nuclei. Even if a neutron is not Bragg reflected, its velocity will generally differ 
from that in empty space and will be approximately given by v = p/m*, where p is the 
neutron momentum and m* (which depends on p) is a neutron effective mass caused 
by the quantum mechanical interaction of the neutrons with the periodic lattice. For 
the momenta satisfying Bragg's law, the effective mass becomes infinitely large to the 
effect that v = 0. Due to the Pauli exclusion principle, free neutrons have different mo- 
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menta p and are therefore diffracted differently. We have systematically calculated the 
momentum-averaged n eutron effect ive mass in all regions of the inner crust using the 
band theory of solids (lChamelll2012r) . The largest effects of Bragg diffraction are found 
in the intermediate parts of the inner crust at baryon densities n ~ 0.02 - 0.03 fm 
where m* ~ 14m„. These results imply that only a very small fraction of neutrons are 
not entrained by clusters and can therefore be considered as effectively "free". Alter- 
natively, these results imply that clusters are effectively much heavier than generally 
assumed (Pet hick et al.ll2010l) . 



4. Conclusions 

The pairing mechanism giving rise to superfluidity in neutron-star crusts is a highly 
non-local phenomenon involving both neutrons bound inside clusters and free neutrons. 
Using the band theory of solids, the critical temperature for superfluidity is found to be 
reduced by ~ 10 - 20% as compared to pure neutron matter. Non-local effects have a 
much stronger impact on superfluid dynamics. 

Due to non-dissipative entrainment effects, some regions of the crust are found to 
strongly resist the neutron superfluid flow. This may have implications for various astro- 
physical phenomena. For instance, the amount of angular momentum that the neutron 
superfluid can possibly transfer to the crus t is limited by entrainme nt, thus challenging 



the interpretation of large pulsar glitches (Cha meT& Carter! 120061) . Because the neu 



tron superfluid is coupled to the crust, entra inment leaves its imprint on the spectrum 
of collective excitations and their damping (IPethick et al.ll201(J : ICirigliano et al.ll2011r) . 
This could affect the thermal relaxation of the crust which has been recently monitored 
in several persistent soft X-ray transients. 
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